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Background of the Invention 
Natural gas Is produced from underground reoorvolr:i 
through wells either alone (gdc wells) or Jn conjunction with 
crude oil (oil wells). Alaoat always associated with gns pro- 
duced from these sources are recoverable liquid hydrocarbons, 
water, and minor amounts of non- hydrocarbon gaseous constituents 
(such as hydrogen sulfide, nitrogen, and carbon dioxide). The 
gas stream Is conditioned and processed In many ways between 
tne time It leaves the well and the various constituents of 
the stream are sold, producing revenue. The result Is that 
the volumes of products sold cannot be directly related to the 
volumes produced from the well without a considerable amount 
of mathematical manipulation of the available volume and price 
data . 

The possible variations from one physical gathering 
system to another are practically infinite in number and scope. 
Added to this is the fact that the gas and products are processed 
and sold under contracts which vary as to their conditions and 
requirements on the seller, the purchaser, anil the processor. 
The combination of these two factors tends to make each alloca- 
tion situation approach unique. The allocation process must 
take into account each and every physical occurrence in the gas 
stream between the producing well and sales point, as well as 
the conditions and requirements of the contracts involved, and 
calculate the equitable share of each product for each entity 
furnishing a part of the total gas stream. Allocation proce- 
dures are dynamic by nature and standardized calculation proce- 
dures, to be useful, must be flexible enough to allow fcr 
changes which take place In the situation, e.g. the addition 
or removal of wells or changes In the gas processing method. 

Prior to the advent of the electronic computer, these 
allocation procedures were performed manually by clerks using 
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calculators and the very tedious nml laborloua calculation 
formulas developed for each Individual situation. In large 
fields the allocation procedure could require in excess of 
40-50 man-hours per month to perform. 

The variability in procedures and the dynamic nature 
of the allocation situation made early attempts at computeriza- 
tion limited successes at best. Each field system was generally 
approached on a single allocation system basis. Even then, the 
dynamics also caused frequent logic changes in the computer 

10 programs to keep up with the physical changes occurring In the 
gas system resulting In programming effort disproportionate to 
the advantages derived. 

The primary objective of the invention is to achieve 
automation through a general computer system which would offer 
the advantages of speed, accuracy, and economy. The results of 
the allocation procedures are available sooner with the cor.puter 
system. The use of the digital computer to make the calculation 
results in the ultimate in accuracy of computation, and the 
desired results are available at lower unit cost than from a 

20 manual system. A further object of the invention is to perr.lt 
universal handling of multiple physical allocation situations 
and dynamics with a general purpose computer system requiring 
no internal programming modification to reflect changes which 
take place in the physical gas gathering and processing complex. 

Summary of the Invention 
This Invention is a system for the allocation of gas 
and liquid products through field gathering and processing. Tnc 
gas stream is carried successively from the producing well 
through various measuring and processing points, all encountered 

30 use or disposition points, and down to the point of entry to a 
gas processing plant or to the point of sale If the gas is not 
processed in a plant. The invention takes as input data the 
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various measured volumes, stream analyses, fc.ctor3, known 
estimates, and formulas for the calculation of vnrloun gar* 
dispositions which take place In the nyntom. Through tho ur.r 
of a defined and coded procedure applying to that particular 
system, It performs all the mathematical calculations and 
manipulations of data required. The output results are the 
allocated volumes of gas and products due to each entity pro- 
viding a portion of the total stream flowing through the field. 
The generalized nature of the system permits Its use in per- 
10 forming these functions for all fields, regardless of the widely 
varied physical makeup and required functions from field to 
field. This is accomplished by the use of a set of generalized 
functions, each of which is coded Into the calculation program 
of the system as an independent section or subroutine, callable 
and performable In unlimited frequency and sequence. These 
general functions cover the spectrum of all types of operations 
necessary in the process of allocation. The set of functions 
which is required to perform any given allocation is coded as 
a special type record of the system master file. This permits 
20 its updating to account for the changes which occur in the field 
operation through the relatively simple procedure of file 
maintenance rather than requiring internal logic changes in 
the computer programs. The actual allocation process is per- 
formed in the calculation program of the system through the 
insertion of input data into a large two-dimensional array 
constructed in the computer core. The function subroutines 
then direct the various calculations, movement and manipulation 
of data in the array until values representing final allocated 
volumes are obtained. The data is then stored on magnetic tapo 
30 or disk devices and reports are prepared through standard data 
processing procedures. Other ordinary data iyrocessing programs 
may be made a part of the overall system to perform functions 
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such as the input of monthly transaction data to the mli'iiu- 
tion program, sorting and printing of various data at several 
points in the system, and normal tests and validation of data. 
Brief Description of the Drawings 
FTO. 1 presents a hypothetical field gas gathering 
system showing some typical sources and dispositions of gas 
which bear directly on allocation procedures. 

PIOS. 2, 3* 5* 6, and 7 present the design logic 
of typical generalized allocation function subroutines used in 
10 the system shown in FIO. 1. 

Description of the Preferred Embodiment (n) 
Raw gas streams pass from oil and/or gas welln Into 
and through gathering systems of diverse constitution and 
orientation. There are numerous possibilities for the trans- 
formation of the character of the stream and the partial use 
or disposition of the stream as it flows through a given 
system. For example, objectionable components may be removed 
from the gas, a portion of the stream may be required for fuel 
to operate various units of the system, and various useful 
20 products may be recovered from the gas for sale. 

The ultimate result is that after all the necessary 
and desirable operations take place, the valuable products of 
the stream are sold or otherwise disposed of, resulting in 
revenue attributable to the producing entities which have 
furnished the gas to the system. The purpose of allocation 
procedures is to assure that the volumes of products and/or 
their value is equitably distributed to the supplying entities 
or, in other words, that each involved owner interest receives 
his fair share of the proceeds. 
30 The invention is described in conjunction with the 

embodiment presented in FIG. 1 which shows a typical field 
gathering system. Gas enters from groups of wells A, B, C, D, 
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E, and P. Note that groups A-D are gas wells while groups E 
and F produce both crude oil and gas. Wells 1-5 are in group A; 
wells 6-9, In group B; wells 10 and 11, in group C; wells 12-18, 
in group D; wells 19-23, in group E; and wells 24-26, in group F. 
The various wells are placed in their respective groups primar- 
ily on the basis of stream pressure. Each group or sub-system 
has to be adjusted to the lowest pressure well present In the 
group. The gas production from each well is measured by flow 
meters m prior to being combined with gas from other wellr. In 
10 the group or with other groups. Flow meters m also measure 

any input or output of gas into or from the system at the other 
locations indicated in FIG. 1. 

Referring to group A, wells 1-5 produce gas at 1000 psi 
that is combined to form stream 29 leading to dehydrator unit }0. 
Fuel required by dehydrator 30 in its operation is taken from 
output stream 31 as indicated by return line 32. Meter m in 
line 31 measures the fuel consumed by dehydrator 30. 

The group B wells produce gas at 500 psi and their 
pressure must be boosted prior to combining with the gas from 
20 the group A wells. The output from wells 6-9 is combined to 
form stream 33 leading to compressor 34 where the pressure Is 
raised to the desired line value of 1000 psi. Output 3tream 35 
from compressor 34 is measured by meter ra and goes to de- 
hydrator 30 where it is combined with the gas from group A. 
Fuel for compressor 34 is taXen from output stream 35 as indi- 
cated by return line 36. Meter m in line 36 measures the fuel 
utilized by compressor 34. 

Following output stream 31 from dehydrator 30, a 
Junction is reached where gas is removed via line 37 for opera- 
10 tion of a drilling rig. Meter m in line 37 measures the volume 
of gas provided the rig. 

The gas continues to move through line 31 to the 
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connection points with lines 38 and 39. Here la a situation 
where some of the owners of this system have agreed to exchange 
gas .with some other system under conditions mutually beneficial 
to both systems. Gas transferred to the other system Is 
measured by meter m on line 38 and gas received is measured 
by meter m on line 39. 

Next, gas produced by the group C wells in gathered 
by line 40 and connects with stream 31. Gas produced by welln 
10 and 11 is already at 1000 psi so no pressure adjustment Is 
10 required • Also, dehydration treatment is not necessary. 

Referring now to group D, wells 12-18 produce high 
pressure, 2000 psi, gas that is subject to a different field 
processing arrangement. Here the gas is combined to form Input 
stream 45 going to low temperature separator 46 which reduces 
the pressure of the gas. Energy thus released la used to 
refrigerate the gas to accomplish dehydration and to condense 
a liquid hydrocarbon product. The condensed hydrocarbon liquids 
pass via line 47 to storage tank 48. Some low pressure, 100 psi, 
gas results as a byproduct of stabilizing the liquid hydrocarbons 
20 and is output through line 49 to line 60 where It goes to the 

second stage of three-stage compressor 55. The main product Is 
1000 psi gas which forms output stream 50. Fuel required to 
operate separator 46 is taken from stream 50 by return lino 51. 
Meter m in line 51 measures the gas utilized by separator 46 
and meter ra in line 50 measures the net volume of 1000 psi gas 
leaving separator 46. Downstream, it will be seen that line 5? 
branches off from line 50. Line 52 carried off gas for use as 
lift gas to aid In oil recovery. 

The remaining segment of the gathering system consists 
30 of gas which has been produced from oil wells. Gas from groups 
E and F will normally be at a considerably lower pressure and 
will contain higher quantities of recoverable liquid hydrocarbons 
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30 



than gas produced by wells A-D. The handling cf the gas and 
its recovery Is generally appreciably different as a result of 
these factors. 

First, it will be noted that gas is introduced into 
wells 19 and 20 via lines 56 and 57 for lifting oil from the 
underground reservoir. This gas is measured by metern m in 
lines 56 and 57 prior to injection into each well. The output 
product from each of wells 19-23 (Including injection gas) is 
combined to form stream 58 going to two-stage separator 59. 
One part of the gas recovered from separator 59 is 100 psi and 
is carried by line 60 to the second stage of compressor 55. The 
remainder of the gas recovered from separator 59 is ?50 psi and 
is carried by line 61 to the third stage of compressor 5*3. 
Meters m in lines 60 and 6l measure the gas volumes thus 
distributed. Recovered liquid petroleum products are carried 
off by line 62 to storage tank 63. 

Production from the group F wells is passed via line 
64 to two-stage separator 65. One part of the gas product from 
separator 65 is 100 psi and is carried by lines 66 and 60 to 
the second stage of compressor 55. The other part of the gas 
product is 250 psi and is carried by lines 67 and 6l to the 
third stage of compressor 55. The liquid residue passes via 
line 68 to oil treater 69 where some 50 psi gas is obtained 
which is transported to the first stage of compressor 55 via 
line 72. The volume of the 50 psi gas is measured by meter m 
in line 72. The residual oil is carried via line 70 to storage 
tank 71. 

After the gas supplied to compressor 55 is compressed 
to 1000 psi, it is carried via line 75 to line 50. Fuel used 
by compressor 55 is indicated by return line 76. Meters m in 
lines 75 and 76 measure the volume of gas introduced into line 
50 and used by compressor 55, respectively. Additionally, some 
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low pressure gas is flared at compressor 55 as shown by metered 
output line 77. 

Finally, line 50 connects to line 31 so that gnn from 
all the wells is Joined together. At valve 41 line 31 divides 
into line 31a goinfc to a gas transmission line and line 31b 
going to a gas processing plant. The entire stream would go 
into one or the other of these branches with meters m measuring 
the volume of gas carried. 

Now let us apply the processes of the present invention 

10 to perform the allocation procedures required in the field system 
shown in FIG. 1. 

First, referring to the group B wells, the total 
volume of fuel consumed by compressor 34 must be allocated to 
each of wells 6-9 since they have gained the benefit of the 
compression. This is done by using the computer coding scheme 
referred to as the "Allocation Function" illustrated in TABLE I. 
By inserting the necessary data into the coding structure for 
the Allocation Function, the total volume of fuel consumed by 
compressor 34 will be charged to each of wells 6-9 in proportion 

20 to the initial volume of gas supplied by each well to compressor 
34. Computer logic for carrying out the Allocation Function 
subroutine is given in FIG. 2. 
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TAMLK I 
ALLOCATION WNCTTON 



MASTER FILE RECORD CODING SCHEME 



ALOC 


K 


L 


C A 


R I C 
K A j °R 
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ALOC = Operating Code 

K =» Number of Allocation Points Involved 

<KX>) 

L » Number of Exception Points to Alloca- 

tion (LSO) 

C A » Column Name Location of Total Volume 

to be Allocated 

R. = Row (Alloc, Pt.) Name Location of 

Total Volume to be Allocated 

C R «= Column Name to Receive Computed 

(Allocated) Volumes 



= Column Name Location of Base or 
Theoretical Volumes 



E,...,E L = Row (Alloc. Pt.) Names of the 
Exception Points 

A 1 ....A„ = Row (Alloc, Pt.) Names of Allocation 
1 * Points 



Referring to FIG. 2, the function and procedures 
performed by the Allocation Function are as follows: 

Input data i3 introduced into blcck 80 to initiate 
the subroutine. The actual volume to be allocated is determine 
in block 8l. If the actual volume is zero, the subroutine is 
terminated by following the YES line from decision point 32 to 
block 9*. If there is an actual volume, the NO line is 
followed to decision point 83. If L=0, i.e. there are no 
30 exception points, the YES line is followed to block 8k. If 

there are exception points to consider, the NO line is followed 
to block 85 where the exception volume is totalized. 

V0L (C RJ E,) = V0L (C R , E,) + V0L (C V , E,) etc ' for (V"- E L> 
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The actual volume (net) Is calculated in block 86 by 

subtracting the total exception volume from the total allocation 
volume. 

ACT VOL » VOL, c R , - [VOL,,, VOL,,, J 

The indexes are changed and the Mar.c Volume accumu- 
lator set to 0 in block 84. The liaac Volume In computed by 
looping through block 87 until all K points have been used. 
BAS VOL - VOL (Cv> A ^ + VOL (Cv> ^ + . . . . VOL^ ^ 

10 When K=0, the YES line is followed from decision point 

88 to decision point 89. If the BAS VOL is something other than 
zero, an allocation factor is computer according to block 90. 

FACTOR » ACT VOL 
FACTOR » HAS VOL 

Three steps take place in block 91. A theoretical 
volume is calculated; this volume is added to whatever voHir.o 
may already be present at the receiving point In the array am! 
the theoretical volume is deducted from the actual volume. 
V0L (C R , A,) - V <*(C R) Al ) + t FACTOR x VOL^ etc. for 

(A^..«.Ag) and 
20 ACT VOL = ACT VOL - T. VOL 

When K=0, the YES line is followed from decision point 
92 to block 93 where the final allocated volume is stored at Its 
code point in the array and the subroutine is terminated. 

An example of the coding required to accomplish this 

follows: 

AIX)C/000V0000/SYSF/1F01/SYSF/MV^^ 
where "1F01" is the code point for the fuel used by compressor 
34 and n 7W06," "7W07, " "^08," and "7W09" are code points for 
wells 6-9. 

30 Thl3 coded statement tells the calculation program 

to call program section ALOC and perform it wiO the inserted 
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data. There are four points for which the allocation la to be 
made and there are zero exception points. The volume to bo 
allocated will be found In the array at the column market! ;;y:*i*' 
(system fuel) at code point 1F0J . The calculated allocated 
volumes are to be placed In the array under SYSF at each or 
points 7V06, 7V07, 7W08, and 7W09. The allocation will be made 
based on the volumes present under the column called KVOL 
(metered volume) for each of the four points Involved. The 
program then runs through the logic for the subroutine deter- 
10 mining the share of compressor fuel chargeable to eacn of well.. 
6-9 and stores thl3 data In the array for future use. 

Turning to dehydrator unit 30 f the fuel metered on 
line 32 must be allocated among all the wells benefiting from 
this treatment. Thus, fuel should be charged to eacn of wells 
1-9 based on the portion of gas they provided to dehydrator 30. 
This would Involve a second charge to the group B wells which 
can be kept separate or combined with the charge relating to 
compressor 34. The same allocation function given above would 
be used except that the fuel would now be divided among nine 
20 points rather than four. The coding for this would be: 
ALOC/0009/0000/SYSF/1F02/SYSF/MVOL/TW01 

/7V02/7V03/AfC^/7V05/^ 
where "1F02" Is the code point for the fuel used by dehydrator 
30 and °7rf01-7WO9" are code points for wells 1-9 . 

Now, referring to the gas output via line 37 to the 
drilling rig, assume that contractual arrangements have been 
made with the owners of well 2 to obtain gas as required to 
operate the rig. The measured volume of gas passed through line 
37 then must be charged to well 2, but not to any other wells 
jO In groups A and B. 

This Is done by use of a coding function called "roint 
Move" which moves an entry from one array point to anot t array 
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point. The Point Move Function la given in TAHLE II and i.hc 
computer logic for it is set forth in FTO. 3. 

TABLE II 
POINT MOVE FUNCTION 



MASTER FILE RECORD CODING SCHEME 



10 



20 



FMOV 


S 


C V 




C R 





PMOV » Operating Code 



Function Sign 



"PLUS" = Value to be 

Added to Value 
in Fid. 
"SUBT" => Value to be 

Subtracted from 
Value in Fid. 
"MOVE " » Value to Overlay 
Value in Fid. 



Cy =x Column Name Location of Data to be Moved 

R v * Row (Alloc. Pt.) Name Location of Data 
to be Moved 

C R =» Column Name Location of Receiving r'ield 

R R « Row (Alloc. Pt.) Name Location of 
A Receiving Field 



30 



Referring to FIG. 3, the function and procedures 
performed by PMOV are as follows: 

Values are entered and indexed according to block ICO. 
The function sign S is then tested at three-way decision point 
101, If S = MOVE, calculations are performed as shown in blcrk 
102 and the volume entry is simply moved to another point in 
the array. If S * SUBT, calculations are performed nr. r.howr* In 
block 103 and the volume entry is moved to another point ~*ntl 
subtracted from the volume already there. If S = PLUS, calcula- 
tions are performed according to block 10k and the volume entry 
is moved to another point and added to the volume already show 
for that point; Whichever route is followed, the subroutine is 
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terminated as shown in block JO'j nftor thn <Jntn nhin. I:: 
accomplished and t.he renu3tc art: :it.onui. 

Following the coding :;ehomo ::hown In TAMLK 11: 
i wy/move/mvol/i r>OJ /:;ALh:/7w<>: , 
where "1S01 M is the code point for the sale of gas via line W 
to the drilling rig and "7V02" la the code point for well 2. 

This coding tells the calculation program to call the 
PMOV section and perform it using the coded data. The volume 
of gas measured on line 37 is entered under "MVOL." The 
computer then transfers this entry to code point "7W02" where 
it is placed under the column called "SALE. ,r The "MOVE" means 
that the entry is transferred, replacing any value already 
present under "SALE." 

Referring next to the gac exenange represented by 
output line 38 and input line 39i assume that the owner of 
well 9 has contracted for and participates In the exchange but 
the other well owners do not. Therefore, any gain or loss of 
gas due to the exchange must be added to or subtracted from the 
output for well 9. 

The Point Move Function Is again used in this situa- 
tion since a simple shift from one array point to another is 
all that is required. Using TABLE II, two separate codings 
are required: 

PM0V/SUBT/MV0L/1E01/EXCH^W09 
and 

PM0V/PLUS/MV0L/1E02/EXCH/7W09 
where "lEOl" is the code point for line 38, "1E02" is the code 
point for line 39, and "TWOg" is the code point for well 9. 

The former coding represents the situation where gas 
leaves the system via line 38 and the latter coding covers the 
case where gau is introduced via line 39. In both cases, 
entries made under "MVOL" are transferred to the code point for 
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well 9* i.e. 7W09, and placed under the column called "EXCH." 
When gas Is being transferred away from the system, the volume 
la subtracted from the value already In the column. When g/is 
Is being transferred to the syctem, the volume la ndJed to the 
volume, If any, already there. 

The wells of group C flow directly to salea by ltne 
31a or to the gas processing plant by line 31B. As a result, 
wells 10 and 11 enter only Into the final allocation procedure 
at the sales point or plant Inlet. (Their effect will be 
10 discussed more fully later on.) 

Referring to the group D wells, fuel required by low 
temperature separator k6 must be allocated to wells 12-18 in 
accordance with the volume Input made by each well. This would 
be done by means of the Allocation Function previously Illus- 
trated using the following coding: 

ALOC/O007/O000/SYSF/1F03/SYSK/MVOL 
/7W12/7W13/7W1U/7W15/7W16/7W17/7W18 
where "1F03" Is the code point for the fuel used by separator X 
and 'TWlS-TWlS" are code points for wells 12-18. 
20 Additionally, the liquid condensate recovered In 

storage tank 48 must be allocated to the wells providing the 
gas from which It was recovered. Since gas generally varies in 
content with respect to those components which make up the ccn- 
densatlon product, It would not be equitable to allocate on the 
basis of gas volume alone. To be accurate one must consider 
the analysis of the gas and calculate a theoretical liquid 
volume based on the analysis with respect to each well. 

The procedure to allocate the liquid volume would first 
Involve the calculation of a theoretical liquid volume from each 
30 well using the "Theoretical Volume Function" represented in 
TABLE III and FIG. k. 
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TABLE III 

THEORETICAL VOLUME CALCULATION FUWCTTfiM 



MASTER FILE RECORD CODING SCHEME 

[ c c| S 



TVOL I K I C„ 



C R 



A«^ » • • . 



TVOL « Operating Code 

K » Number of Allocation Points 

Involved (K>0) 

C c » Column Name; Location of 

Factor to be Used in Calc. 

C v = Column Name; Location of Data 

to be Used in Calc. 

C R m Column Name to Receive the 

Calculated Values 

A 1 ....A JC e Row Names of the Allocation 
Points Involved 



Referring to FIG. 4, values are tabled and indexed in 
block 108. At decision point 109, determination is made re- 
garding the value of K. If K«0, then the subroutine terminates. 
If Wl, then calculations are performed according to block 110. 
The theoretical volume is calculated by multiplying the original 
data by a factor and adding the result to any value already 
there. 

V0L (C R , A X ) « V0L (Crj MV0L (v v x FACT0R (Ccj ^ 

The above step is then repeated for points A 2 ....A K as indicated 
by block 111. 

The coding required to perform the Theoretical Volume 
Function in regard to the condensates recovered in tank 48 is 
as follows: 

TVOL/0007/SEPT/^VOL/CONP/7W1 2 
/7Wl3/7wii;/7wi5/7Wl6/7Wi7/7Wl8 
where "SEPT" is the array code for the separator GPM factor and 
"7*12-7*16" are code points for the allocation locations. 
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This coding calls the program section called TVOL and 
has the computer perform it using the coded data. Volume 
calculations are to be performed for five different points 
using the analysis factors given under the array point SEPT. 
The theoretical condensate production (CONP) io then calculated 
by multiplying the metered gas volume for each well (MVOL) by 
its analysis factor* 

Following the performance of the theoretical volume 
function as set forth above, the actual recovered condensate 
10 liquid would be allocated to the wells of group D using the 
previously illustrated Allocation Function. To complicate 
matters, assume that wells 17 and 18 have contracted to receive 
their full theoretical volume of condensate. The coding would 
then be: 

AL0C/0007/0002/C0KP/1X01/TMP1/C0NP 

/7V17/7V18/7W12/7W13/7W14/7W15/7W16 
where "1X01" is the point name of the actual stored condensate 
and "TMP1" is a temporary storage column. 

Wells 17 and 18 are listed as exception points and 

20 will automatically be eliminated from the allocation procedure, 
with the .full theoretical volume replacing the allocated volume 
which would have been calculated. 

The low pressure gas which Is fed to compressor 55 by 
line kg must al3o be credited to the wells which furnished the 
liquid product. Therefore, it will be allocated on the came 
basis as the liquid condensate, i.e. the allocation function In 
used based on the theoretical volume calculations. The coding 
would be the same as that given above, except for the column 
location of the gas to be allocated. 

30 Downstream of separator 46, lift gas is removed via 

line 52. This volume of gas must be allocated to the well? 
that feed separator 46, I.e. wells 1P-16. Allocation in thin 

-17- 



Instance will be based on the gas stream going lo separator 
rather than the theoretical liquid volume, using the Allocation 
Function. Coding would be: 

AU)C/0007/0000/^OL/1L01/GLFT/MVOL 
/7W12/7W13/7W1V7W15/7W16/Trfl7/7W18 
where "lLOl" is the code point for the lift gas witndrawn 
through line 52. 

This coding tells the computer to perform f, ALOC M 
using the entry made at code point 3L0] under the column cnlled 
MVOL (for gas lift). The original volume entry Is to be 
allocated among code points 7W12-7W18 and stored under the 
column called GLFT using the metered volumes MVOL as the basis. 

Referring to the group E wells, gas utilized by wells 
19 and 20 must be charged to wells 19 and 20 and Kept separate 
from the gas produced by the other group E wells. This Is 
accomplished by a function called "Gas Lift" which is illus- 
trated in TABLE IV and FIG. 5. 
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TABLE IV 
OAS LIFT FUNCTION 



MASTER FILE RECORD CODING SCHEME 



GLFT 



K| C £ C L |C. 



OUT 
K • 



- Operating Code 



>^ • • • 



• • • 



» Number of Allocation Points 
Involved (K>0) 

# There must bo a rorrn spend !ng 
point In the "r" net for ea^n 
point in the W A" -oi. 

= Column Name; Location of Volume 
Entering Well 

= Column Name; Location of Volure 
Leaving Well 

= Column Name; Receiving Location 
When Difference Is Positive 

= Column Name; Receiving Location 
When Difference Is Negative 

• A K = Row Names of Allocation Points 
Involved 

= Row Names of Gas Lift Measuring 
Points Corresponding to the "A 
Set 



Referring to FIG. 5, the operational sequence 
initiated by entering and indexing the values snown in z 
115. If K=0, the YES line is followed frc= decision pci- 
and operations are terminated. If :<*e, tnen tne in terror 
volume (called "IVOL") is calculated in block 117 accord 
the following equation: 

IVOL = VOL (Cl> v - VCL (Cs> 

If IVOL is positive, the *r£S line is followed - 
decision point 118 to block 120 where tne following cal:. 
are made: 



V ° L (C + , Bl ) = V0L (C^ B.l * *VCL 
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If IVOL is zero or negative, the NO line is followed 
to block 119 where the following calculations are made: 
V0L «L, B x ) - V0L ( C _, Bl ) + IV0L 

The outputs from blocks 119 And IPO both lead to 
block 121 where indexing is carried out 30 that, the next, 
allocation point will be considered, i.e. the procedure in 
repeated for points Ag....A K . 

An example of the coding would be: 
OICT/0002/OTOI^OiyFROT^ 
10 where "7019" and 7020" are code points representing gas injec- 
tion into wells 19 and 20, respectively, and 7W19 and 7W20 are 
code points for wells 19 and 20. 

This coding tells the calculation program to call the 
subroutine OLFT and perform its logic with the above data. 
This in effect compares the volume of gas entering with the 
gas recovered from each well. Tf the latter la greater than 
the former, the difference is placed In the column called FRMT 
(formation gas). If the opposite Is true, the negative amount 
is placed in the column called NJCT (injection gas). 
20 The gas from the group E wells Is recovered by means 

of separator 59 and is allocated to wells 19-23 on the basis 
of their gas-oil ratios. These factors are used with the 
Theoretical Volume Function previously discussed to determine 
the theoretical volurue of gas produced by each well. This 
volume then becomes the theoretical basis for allocating the 
gas recovered by separator 59 in streams Go and 6l using the 
Allocation Function. 

Referring to the group F wells, gas obtained from 
separator 65 in streams 66 and 67 is allocated to wells <?4-26 
30 on a theoretical volume basis. First, the Theoretical Volume 
Function is used based on the gas- oil ratio of each well; then 
the Allocation Function is applied. 

-P0- 
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A special consideration in this segment of the system 
would be the gas (vapors) obtained via line 72 from oil treater 
69. This gas is allocated on the basis of oil volume from eacii 
well involved since the oil is the source of the gas. The 
standard Allocation Function is used. 

The gas from the group E and F wells, after leaving 
separators 59 and 64 and oil treater 69, flows to compressor 55 
which boosts the pressure to the point where the gas may enter 
the reat of the system. Gas is also supplied to compressor 55 
10 from the group D wells via stream 49. Here, again, fuel is 

utilized in the compression operation which must be allocated 
and charged against the wells furnishing the gas, i.e. wells 
12-16 and 19-26. The theoretical volume calculations previously 
referred to are used with the Theoretical Volume Function to 
determine the theoretical volumes of gas produced by each well. 
This then becomes the theoretical basis for allocating the fuel 
used by compressor 55 using the Allocation Function. Gas flared 
at outlet 77 is similarly allocated and charged to the con- 
tributing wells. 

20 At point all the Inputs and dispositions of gas 

indicated on FIG. 1 have been considered except for two final 
operations. Let us assume that the gas Is sold directly to a 
transmission line. The gas throughput measured by meter n In 
line 31a then represents the volume actually sold to the 
purchaser. This volume may not agree with the mathematically 
determined total after all of the system operations are 
completed so a final allocation must be made to determine earn 
well's share of the sale volume. This is accompllsned by use 
of the "Allocation Base Function" illustrated In TABLE V and 

30 FIG. 6. 



OliKWENT PU# T CATIOrsi 



mm 



9591 69 



10 



TABLE V 
ALLOCATION BASE FOHCTIOW 

RASTER PILE RECORD CODING SCHEME 



ABSE 


K 




A^ • • • • Aj^ 



ABSE » Operating Code 

K » Number of Allocation 

Points Involved (K>0) 

C v = Column Name; Location 

of Peglnnlng (Cross) 
Volume 

A«....A„ - Row Names; Allocation 
Points Involved 



20 



Referring to FIG . 6, values for K, C y , and A are 
entered as shown In block 125. If K^C, trie YES line Is followed 
from decision point 126 and the subroutine Is terminated. If K 
Is something other than zero, the NO line Is followed to bloc* 
127 where the Initial Index column Is called. In block 128, 
the beginning volume Is placed In the receiving position for 
point A and then the Initial column Index Is compared wltn the 
final column Index. If they are not equal, the following 
calculation Is performed as per block 131. 

V0L (BASE, Al ) " VCL (BASE, A,) * V0L ( Xl . A.J * V0L (X„. A,) 



where X, .... X 
l n 



n' 1 

the disposition ir.d receipt colunis of the 



30 



array. 

When the Initial Index becomes equal to the final 
Index (block 132), the YES branch Is followed from the decir-ion 
block which advances the Index counters and repeats the rrocer.;. 
for points A£-**-Ag. 

The coding for the Allocation Base Function would be 
as follows: 



UJbKWJtlNT PU15Ll#tT10N&' 



9591 69 

ABSE/0026/W0L/7WOl/7W0?/7WO3/7rf0V^ 

/7W08/rrf09/mo/mi/7W12/r^l3/^ 

/7wi8/7wi9/7W20/rrf2i/rrf2?/^ 
This coding tells the calculation program to call the 
Allocation Base Function and perform It with the above dntn. 
The program would take for each well in the system itr. metered 
gas volume and deduct from this each of the disposition;-, of gns 
chargeable to that well and come up with a remainder volume 
which represents theoretically what that well has available to 
sell at line 31a. These volumes then constitute the base on 
which the final allocation of sales volume is made to each 
contributing well using the actual metered sales volume in the 
Allocation Function. 

Once the allocated sales volume for each well is 
determined, it may be desirable to determine each well owner's 
share of that gas* The "Ownership M Function shown in TABLE VI 
and FIG. 7 can be used for this purpose. 



TABLE VI 

OWNER INTEREST CALCULATION FUNCTION 
MASTER FILE RECORD CODING SCHEME 



OWNR 


K 


c v 


A^. . . «Aj^ 



OWNR * Operating Code 

K = Number of Allocation 

Points Involved (K>0) 

C„ - Column Name; Location 

of Volume to be Split 
Among Owners 

A,....A K = Row Names; Allocation 
Points Involved 

Referring to FIG . 7, the values for K, C y> and A are 
tabled in block 135. The subroutine terminates at decision 
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point 136 if K=0. If K Is not equal to zero, then some 
additional initializing of values and Indexes takes place in 
block 137. The next decision point checks to see if the owner 
factor is zero. If it is, block 139 increases the index *by one 
and following a test to assure that the indexed owner factor 
is not greater than the final owner factor (block l4o), the 
process is repeated. If the condition of 140 tests true, 
indexes are increased and the routine returns to beginning, 
repeating for point etc. 
10 In the case where the tested owner factor of block 138 

is something other than zero, the logic proceeds to a further 
test (block 142) of the indexed owner factor compared to the 
final. If the test is passed at this point, the Indication is 
that there is only one owner and the final calculation of block 
143 is made and the indexes advanced to return for points A 2 , 
etc. 

If the test of block lk2 fails, then further proc- 
essing is indicated. The index of the owner factor is moved 
to a previously used index in block 144 and increased by one 
20 in block 145. 

Then the owner index test Is made again (block 146) 
with the same results as previously if the test passes. 

If the test falls, this means that there are other 

owners to be considered and the calculation is made according 

to block 147. An intermediate volume is calculated by the 

following formula: 

ALOC VOL X OWNER FACTOR , . = I VOL, 

n+1 1 

This calculated volume is then added to the owner volume already 
present for owner n+1 and a second intermediate volume is calcu- 
30 lated which is an accumulation of the I VOL^ calculated above. 

This is used in the final calculation of block 143 to determine 
the first owner's share by difference. 

-24- 
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Indexing Is then performed as ?hown in block ]4l go 
that the above sequence is repeated for points A 2 ..,.A K . 

The Ownership Function would be coded as follows: 

0VW0026/AV0iy7W01/7W02/W 
/7W08/7W09/7W10/7Vll/7Vi2m 
/7Wl8/7W19/TW20/7W21/7W22/7W23/7W24/rrf2 

The participation factors of the various owners In 
each well were previously entered into the master file of the 
computer and are available for recall. In this instance the 
10 coding directs the calculation program to perform the Ownership 
Function with the Indicated data. This would result In the 
division of the final allocated gas volume (AVOL) for each well 
in the system to each separate well owner in accordance with 
the percentage ownership which he h/13 in the well. 

Calculations similar to the above would be performed 
if the gas were going to the gas processing plant Instead of 
directly to a transmission line. The flow through line 31b 
would be measured by meter m and the recorded volume allocated 
to each contributing well using the Allocation Base Function 
20 and the Allocation Function. 

While the aforegoing has illustrated the application 
of the present invention to one specific embodiment, it is 
obvious from a consideration of the design concepts of the 
system that it is open-ended from the standpoint of the incor- 
poration of additional general mat.^ematlcal functions Into U\o 
logic of the system. This may be done at pny time a need ari: n: 
which cannot be handled by the existing functions. It lc also 
obvious that the number of allocation points which may be 
handled is also capable of being increased or decreased in 
30 any number, the only limit being the maximum size of the 
computer storcge available. 

In addition to the Allocation Function, the Point 
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Move Function, the Theoretical Volume Calculation Function, the 
Gas Lift Function, the Allocation Base Function, and Owner 
Interest Calculation Function already discussed, other general 
functions which have been employed by Applicants in regard to 
particular field gas systems Include the following: 

The "Volume Split Function" which divides composite 
volumes into two or more components based on ratios of certain 
special factors. 

The "Column Move Function" which moves data items 
10 from one column to another In the array. 

The "Row Move Function" which move3 data from one row 
to another in the array. 

The "Column Balance Function" which subtracts 
designated point values from a beginning value and stores the 
remainder at a selected point. 

The "Row Balance Function" which subtracts a desig- 
nated column value of a point from a beginning column value 
and places the answer In a selected location. 

The "Allocation of Fuel Function" which is a special 
20 case of the Allocation Function in which the volume of fuel to 
be allocated is first calculated from a set formula. 

The "Row Summation Function" in which the sum of 
various column data items appearing in one row is calculated 
and placed at a designated location. 

The "Column Summation Function" in which the sane 
function is performed for data Items In one column of many rews 
or points. 

The "Pressure Base Conversion Function" in which the 
standard pressure basis of gas measurement may be changed to 
33 conform to variations required by contract, state, or company 
differences . 

The "Zero Function" which permits the elimination of 
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data from designated location:; In Uu» army wli.hln limit.::. 

The "Assignment Function" which permit;! the special 
assignment or allocation of volumes to designated points 
depending upon predetermined conditions. 

The "Multiple Element Row Balance Function " .which 
performs the same function as the Row Balance Function except 
that more than one designated element may be subtracted. 

The "Compare and Assign Function" which compares cne 
value with another and makes calculation and assignment based 
on the result of the comparison. 

The "Partial Owner Calculation Function" which perr.lt 
the calculation of owner shares in an entitles volume on a 
designated partial basis. 

The present invention is also applicable to situation 
other than field gas allocation systems. For example, the 
invention can be used to great advantage in the case of gas 
processing plants. Many things occur within the confines or a 
gas processing plant which change the nature, quality and/or 
quantity of the gas stream being processed. All these events 
can be accounted for in allocation procedures similar to those 
presented for the field gathering system regardless of the 
nature of the plant or the conditions governing the processing 
of the gas and handling of the resulting products. 

The allocation functions previously presented in 
conjunction with the field gas allocation system can be applied 
to a gas plant system, subject to certain modifications in 
array size, data content, and functions performed which will 
be obvious to those skilled in the art. (However, there Is no 
need for the Gas Lift Function or the Allocation of Fuel 
Function.) Additionally, certain general functions which are 
uniquely associated with plant activities are given below: 

The "Value Calculation Function" which permits the 
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calculation of the money value of various volumes of product 
and gas and the storage of these values in selected locations. 

The "Tax Calculation Function 1 ' which performs the 
computation of various forms of taxes which must be deducted 
from gross values prior to determining the actual allocated 
values due leases or allocation points. 

The "Hydrocarbon Shrinkage Calculation Function" which 
calculates from Input factors the amount of nhrlnknge resulting 
in the gas stream from the removal of various hydrocarbon 
products from the stream. 
10 The "Inventory Status Calculation Function" which 

calculates the position of inventory relative to sales and 
production and permits the determination of the source of 
disposed volumes of product and thus the price which is appli- 
cable. 

The "Inventory Balance Function" which then takes the 
status and sales and production and calculates a new closing 
inventory for each product of the plant. 

The present allocation systems have been coded in 
the Cobol computer language for use on IBM 360 computers, but 

20 It is obvious that the invention can be extended to other 
computer languages using any desired hardware and software 
systems. The complete software system consists of several 
conventional data processing programs for producing reports, 
sorting data, and changing record formats. In addition to 
these are the three programs which comprise the heart of the 
system by Incorporating the concepts which have been explained. 
These are the master file maintenance program which keeps the 
data current and permits changes In the allocation procedures 
through file updating procedures. The second is the Indexing 

30 program which provides the means of Indexing or locating all 
data items within the computational array. The third is the 
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computational program itself which performs nil the required 
manipulations and calculations and produces nn output record 
of the results. 
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computational program Itself which performs nil the required 
manipulations and calculations and produces an output record 
of the results. 
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The embodiments of the invention in which an exclu- 
sive property or privilege is claimed are defined as follows: 

1. A system for allocating gas and liquid products 
in conjunction with a field gas gathering and processing system 
comprising: 

(a) a digital computer operable to solve alloca- 
tion functions physically represented by electrical signals 
within said computer, said functions representing uses and dis- 
positions of products encountered in said field gas gathering 
and processing system, and the relationship of the uses and 
dispositions to individual allocation points, 

(b) means for metering gas flow at points of 
allocation and points of use and disposition in the gas gather- 
ing and processing system, and for generating an electrical 
signal applied to said digital computer rolatable to the 
metered gas flow, 

(c) means for generating within said digital 
computer signals representing a two-dimensional array, alloca- 
tion point index terms being defined in one dimension and use 
and disposition index terms being defined in another dimension, 

whereby said digital computer performs the 
mathematical calculations prescribed by the allocation functions 
in respect to the signals relatable to the metered gas flew, and 
stores signals representing results of the calculations in said 
array under appropriate index term.;. 
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